Zeocin-induced DNA damage response in barley and its dependence on ATR

DNA damage response (DDR) is an essential mechanism by which living organisms maintain their genomic stability. In plants, DDR is important also for normal growth and yield. Here, we explored the DDR of a temperate model crop barley (Hordeum vulgare) at the phenotypic, physiological, and transcriptomic levels. By a series of in vitro DNA damage assays using the DNA strand break (DNA-SB) inducing agent zeocin, we showed reduced root growth and expansion of the differentiated zone to the root tip. Genome-wide transcriptional profiling of barley wild-type and plants mutated in DDR signaling kinase ATAXIA TELANGIECTASIA MUTATED AND RAD3-RELATED (hvatr.g) revealed zeocin-dependent, ATR-dependent, and zeocin-dependent/ATR-independent transcriptional responses. Transcriptional changes were scored also using the newly developed catalog of 421 barley DDR genes with the phylogenetically-resolved relationships of barley SUPRESSOR OF GAMMA 1 (SOG1) and SOG1-LIKE (SGL) genes. Zeocin caused up-regulation of specific DDR factors and down-regulation of cell cycle and histone genes, mostly in an ATR-independent manner. The ATR dependency was obvious for some factors associated with DDR during DNA replication and for many genes without an obvious connection to DDR. This provided molecular insight into the response to DNA-SB induction in the large and complex barley genome.


Establishing barley in vitro DNA damage treatment
We developed a protocol to assess barley plant sensitivity to DNA damaging treatments.Approximately 20 dissected mature barley embryos were cultured on 100 ml of regular ½ MS medium in plastic cultivation boxes, resulting in synchronized germination approximately 24 h after initiation of cultivation.
Collectively, we observed reduced root growth and shorter stems after zeocin treatment, with 100 μg/ml zeocin being sufficient for significant changes.Thus, we successfully established zeocin-induced in vitro DNA damage conditions in barley.

Hvatr.g plants are sensitive to zeocin treatment
To explore barley's DNA damage response (DDR) in the context of DNA damage signaling, we utilized a loss-offunction mutant allele of ATR kinase (hvatr.g) in the Sebastian background 21 .Previous study on hydroponically grown hvatr.gplants demonstrated shorter roots compared to cv.Sebastian (WT) 17,21 .We observed a similar trend in plants grown from rescued embryos (n = 9-10) on solid media in vitro (Fig. 2A).Cultivar Sebastian was used as wild-type, and its phenotypic responses were comparable to those of cv.Golden Promise.Under mock conditions, hvatr.gplants exhibited a 60% reduction in root length compared to WT plants (13.4 ± 0.3 mm and 32.8 ± 1.4 mm, respectively) (Fig. 2B).Following 100 μg/ml zeocin treatment, the average root length was 15.8 ± 0.8 mm for WT plants and 8.9 ± 0.1 mm for hvatr.gplants (Fig. 2B).This corresponded to a 51.5% reduction for zeocin-treated WT plants and a 33.8% reduction for zeocin-treated hvatr.gplants (Fig. 2B).The greater reduction in WT plants likely indicates their higher potential for root shortening compared to hvatr.g plants, which already had significantly reduced root length.Similar results were obtained for cumulative root length and the longest root (Supplemental Fig. 1).The median root number was significantly higher in hvatr.gplants (median root number = 11) compared to WT plants (median root number = 8) (Fig. 2C).After zeocin treatment, WT plants maintained the median number of roots (median root number = 8), while the number decreased (median root number = 7; 36% reduction) in hvatr.gplants.Additionally, we analyzed plant height (Fig. 2D).These findings suggest that total shoot length is the most informative characteristic for assessing the DNA damaging effects of zeocin on barley hvatr.gplants.

Identification of barley homologs of DDR genes
To prepare for the transcriptomic analysis, we identified barley homologs of known DNA damage response (DDR) genes.Using the published list of 321 Arabidopsis DDR genes 22 , we performed a series of BLASTs and selection steps, resulting in the discovery of 421 barley homologs (Supplemental Dataset 1).Visual inspection revealed the absence of BRCA1 due to a lower homology.This gene was added manually and made the final set of 421 barley homologs.Approximately 50% of the identified proteins (n = 204) were coded by single-copy homologs, such as ATM, ATR , WEE1, BRCA1, BARD, and RAD51.We also observed an expansion of PCNA2ed family to 11 members in barley.Arabidopsis PCNA2 interacts with the translesion synthesis polymerase POLH and is involved in the repair of UV-induced damage 23 .www.nature.com/scientificreports/Of special interest was the master regulator of plant DDR SOG1 and its potential homologs.SOG1 is a member of a large NAC (NAM, ATAF1/2, and CUC2) family of transcription factors in plants and a list of 167 barley NAC domain containing proteins was published, highlighting HORVU.MOREX.r3.7HG0670800 (alias HORVU7Hr1G042420) as the most likely homolog of Arabidopsis ANAC008 (AtSOG1) 24 .On contrary, another gene HORVU.MOREX.r3.6HG0590960 (alias HORVU6Hr1G053540) was suggested as HvSOG1 based on upregulation upon DNA damage treatment 25 .Therefore, we analyzed SOG1 in barley in more detail.Via BLAST we found five barley genes that shared a significant homology with the AtSOG1 (Supplemental Dataset 1).Three were analyzed previously 24 .To assess their relationships, we performed a phylogenetic study using five SOG1 and SOG1-LIKE (SGL) factors from barley, two from rice, two from maize and three from Arabidopsis (Fig. 3).The Arabidopsis factors included AtSOG1/AtANAC8 and two closely related transcription factors AtANAC44 and AtANAC85 that were also shown to play a role in DNA damage response 26 .
Multiple sequence alignment and phylogenetic tree revealed three main clades, (i) the AtANAC005 and AtANAC006 outgroup, (ii) the SOG1 and (iii) the SGL clades.The SOG1 clade contained Arabidopsis SOG1, rice OsSOG1, both maize SOG1 candidates (ZmNACTF99 and ZmNACTF08) and HORVU.MOREX.r3.7HG0670800, suggesting it as the most likely barley SOG1 ortholog (Fig. 3).The SGL clade, defined based on a lower 45% confidence at its root but high confidence at the branches, contained AtANAC44 and AtANAC85, rice OsSGL, and the four remaining barley candidates.The barley protein HORVU.MOREX.r3.6HG0590960 was most closely related (63%) to the OsSGL; therefore, we named it as HvSGL.A.The remaining three candidates HORVU.MOREX.

Transcriptomic responses of WT barley to zeocin treatments
We analyzed the transcriptional response of barley plants to zeocin treatments.Initially, we examined the expression of barley BRCA1 and RAD51 genes in seven-day-old WT (cv.Golden Promise) plants treated with 300 μg/ ml zeocin for 0, 0.5, 1, and 6 h using RT-qPCR.Both BRCA1 and RAD51 showed significant up-regulation, with BRCA1 exhibiting a six-fold increase and RAD51 showing an almost 11-fold up-regulation at the 6-h time point (Fig. 4).Subsequently, we conducted RNA-seq analysis of root apices isolated from cv.Sebastian WT and hvatr.gplants treated with 0 (mock) and 500 μg/ml zeocin for 6 h, aiming to identify differentially expressed genes (DEGs).
In zeocin-treated versus mock-treated WT plants, we identified 719 DEGs, including 404 significantly up-regulated genes and 315 down-regulated genes (Supplemental Dataset 2).Gene ontology analysis of the up-regulated DEGs revealed their involvement in biological processes such as glutathione metabolic process (GO:0006749), response to stimulus (GO:0050896), response to stress (GO:0006950), and defense response (GO:0006952) (Fig. 5A; Supplemental Dataset 3).Notably, many up-regulated DEGs were associated with oxidative stress signaling (Fig. 5B), including 14 GLUTATHIONE S-TRANSFERASEs (GSTU) and 14 UDP-GLUCOSYLTRANSFERASEs (UGT ) involved in reactive oxygen species detoxification.Other oxidative stress-related genes, such as OXIDA-TIVE STRESS 3, reductases, and a homolog of CATION EXCHANGER 5, were also identified.The observed oxidative stress response can be attributed to the radiomimetic activity of zeocin, which generates hydrogen peroxide upon interaction with a metal ion, leading to oxidative damage to DNA and other cellular components 27 .
We found 7.6 and 2.8-fold up-regulation for RAD51 and BRCA1 in zeocin-treated WT plants, respectively, indicating activation of DDR also in our RNA-seq experiment.Subsequently, we inspected the full lists of significantly down-regulated and up-regulated DDR candidates.In accordance with the phenotypic data, zeocin treatment had a profound inhibitory effect on the cell cycle, including DNA replication and mitotic cell division.We found significant down-regulation for six CYCLINs (CYC ) (CYCB1;2.A-C, CYCB1;3.A-B; CYCB2;3.B) and two CYCLIN DEPENDENT KINANSES (CDKs) CDKB1;2 and CDKB2;2, whose Arabidopsis homologs promote G1/S and G2/M transitions, respectively.
Although the GO term analysis indicated genes related to DNA replication are both up-and down-regulated, a more detailed analysis revealed repressive effects of zeocin treatment on DNA replication.Seven genes marked by the GO term analysis as involved in DNA replication were up-regulated by zeocin treatment: RPA1, RPA2, RAD51, DNA POLYMERASE ZETA SUBUNIT REV3, and RIBONUCLEOSIDE-DIPHOSPHATE REDUCTASE.We argue that the up-regulation of these genes should not be taken as indication for active replication but rather as a sign of its repression and ongoing DDR.The RPA subunits play roles in both DNA replication and repair, DNA POLYMERASE ZETA and RAD51 are primarily associated with DDR, and RIBONUCLEOSIDE-DIPHOSPHATE REDUCTASE is involved in nucleotide synthesis.Next, we found down-regulation of almost all genes coding subunits of the Minichromosome maintenance complex (MCM), DNA primase and REPLICATION PROTEIN A complex.The reduced replication was recognizable also by looking at the expression of histones and replication-coupled chromatin modifiers.In total 63 HISTONE genes were down-regulated, including 15 copies of H2A in several variants, 9 copies of H2B, 23 copies of H3.1 and 16 copies of H4.The downregulated   PROTEIN 6).There were also several DDR-associated down-regulated genes, most notably barley homolog of BARD1.
Zeocin up-regulated genes pointed to specific directions.The up-regulation of two PARP genes and MRE11 homologs indicated presence of both SSBs and DSBs, respectively.We also saw up-regulation of three RPAs: RPA70C.C, RPA2.K and RPA1A that suggest an increased amount of single stranded DNA.At the level of DDR signaling, we found up-regulation of three SGL factors SGL.A, SGL.D, and SGL.C, but not SOG1.The upregulation of barley homolog of Arabidopsis CDK inhibitor SMR3 is in agreement with the alterations in cell cycle and enhanced levels of endoreduplication.The other responses, presumably downstream of SOG1/SGLs

Effects of hvatr.g on gene expression under mock conditions
We examined the effects of ATR mutation on transcriptome by comparing expression in hvatr.g (Fig. 6).The hvatr.g plants showed a significant (adjusted-p ≤ 0.05 and a log 2 fold change ≤ − 1 or log 2 FC ≥ 1) up-regulation of 1150 and down-regulation of 1457 genes relative to WT plants (Supplemental Dataset 2).Gene ontology analyses of the up-regulated genes suggested a connection to only a few categories of biological processes, most related to protein production and/or modifications, as in GO:0006749 (glutathione metabolic process), GO:0006468 (protein phosphorylation), or GO:0006575 (sulfur compound metabolic process).Additional gene categories include xenobiotic export and transport from cell (GO:0046618, GO:0042908).Categories of down-regulated processes show a miss-regulation of transmembrane transport, especially concerning nitrate (GO:0015706, GO:1902025) (Fig. 6A, Supplemental Dataset 3).
Among the most significantly upregulated transcripts was 64-fold upregulated TFIIIA (Fig. 6B), which is known to regulate 5S rDNA transcription.Furthermore, hvatr.gplants had over 182-fold higher expression of The hvatr.g plants showed also up-regulation in genes of the DDR pathway like SGL.C, suggesting increased genome instability even at mock conditions.Furthermore, there was an up-regulation of the MutS homolog 5 (MSH5) that was implicated in mismatch repair in Arabidopsis 28 and nucleotide excision repair and proteolysis associated factors RADIATION SENSITIVE 23 (RAD23) and CELL DIVISION CYCLE (CDC48).The strongly down-regulated genes in hvatr.g were eight copies of NICOTIANAMINE SYNTHASE (Supplemental Dataset 2).These genes should be strongly expressed in roots, where it regulates intake of iron 29 .In rice it was shown to be under direct control of NAC-family proteins during drought stress 30 .The same pathway in rice up-regulates the expression of genes involved in membrane modification genes and transport 30 .Even more remarkable was the decrease in the expression of transmembrane transporters.Among the 71 down-regulated transmembrane transport-related genes, 25 are associated with the transport of nitrogen-based compounds (Supplemental Fig. 3B).Comparison of up-regulated and down-regulated genes with products involved in nitrogen transport or metabolism, confirmed this further (example genes Supplemental Fig. 3A, full gene list Supplemental Dataset 4).The most prevalent down-regulated genes in this group were HIGH AFFINITY NITRATE TRANS-PORTERS and NRT1/PTR FAMILY TRANSPORTERS (6 genes each).
In summary, hvatr.gplants show a well notable pattern of transcriptional changes that includes a mix of DDR, transcriptional and translational responses.The groups of up-and down-regulated genes also indicate enhanced turnover of genic products from the transcription, through protein synthesis, modification and degradation.

Differences in reaction to zeocin caused by hvatr.g mutation
To assess the DNA damage response in plants lacking the functioning ATR kinase, we looked on differentially expressed genes in zeocin and mock-treated hvatr.gplants.There were in total 424 genes up-regulated and 622 genes down-regulated in zeocin-treated compared to mock-treated hvatr.gplants (Fig. 7A, Supplemental www.nature.com/scientificreports/Dataset 2, Supplemental Fig. 2).Gene enrichment analysis showed categories similar to those described in zeocin-treated WT plants (Supplemental Fig. 3).To filter for the ATR-specific responses, we visualized the data as Venn diagrams of genes with significantly changed expression by zeocin-treatment in WT and hvatr.gplants (Fig. 7B, Supplemental Dataset 5).The majority of zeocin-induced transcriptional changes in WT were ATR-independent with 260 genes upregulated (64.4% out of 404) and 264 genes down-regulated (83.8% out of 315) (Fig. 7B, WT and overlap).This included all major up-regulated DDR genes related to SSB repair (PARP), translesion synthesis (REV3), direct hydrolysis (TDP) and HR (MRE11, SGL.D, SGL.C, SMR3, RAD54, BRCA1, BCP1, BCP4, GMI1, NSE4.A).Similarly, genes down-regulated by zeocin treatment in ATR-independent manner included many positive regulators of cell cycle, replication factors and histones (Supplementary Dataset 4).When looking on dependence of transcriptional change on ATR upon zeocin treatment, we found 144 up-regulated genes (35.6% out of 404) and 51 down-regulated genes (16.2% out of 315).The DDR associated genes up-regulated in ATR dependent manner included RECQ HELICASE 3 (RECQL3) and DNA POLYMERASE DELTA catalytic SUBUNIT 1 (POLD1).Both proteins are involved in processing lagging strand during replication where RECQL3 unwinds it and shows a substrate preference for nicked Holliday junctions 31 and POLD performs its synthesis 32 .
Finally, the last category represented zeocin-induced transcriptional changes unique to hvatr.g plants.There were 358 such down-regulated genes among which dominated 67 histone genes of all types (2× H1, 17× H2A in different variants including two copies of H2A.X variant associated with marking the DNA damage sites, 17× x H2B, 5× H3.1 and 26× H4).The down-regulated replication factors included large subunit of DNA primase, MCM3 and GINS complex subunit PSF1, and cell cycle genes MAD2, AURORA KINASE 1.A (AUR1.A) as well as microtubule-associated protein TORTIFOLIA1.From the DDR genes, we found significant down-regulation of RECQL2 helicase, RecQ-mediated genome instability protein 1 (RMI1), two uncharacterized DNA ligases and F BOX-LIKE17 (FBL17).The FBL17 is interesting candidate as its loss of function mutant was identified in an atr phenotype suppressor screen in Arabidopsis and was described as SUPPRESSORS OF ATR 1 (SOAT) 33 .When inspecting 164 genes up-regulated by zeocin treatment in hvatr.g and not in WT, we found APURINIC ENDONUCLEASE-REDOX PROTEIN (ARP) that is a major endonuclease involved in base excision repair (BER), RAD23 that has a role in NER and interacts with the 26S proteasome components, and also ESSENTIAL MEI-OTIC ENDONUCLEASE 1B (EME1B) that is known to interact with endonuclease MUS81.Up-regulation of these components indicates specific shift in the use of DDR pathways towards nucleolysis, BER and NER and in hvatr.gplants.

Discussion
We established in vitro DNA damage treatment conditions for barley, generated a list of barley DDR genes including the SOG1 and SGL family and performed transcriptomic analysis of WT and hvatr.gplants in response to zeocin treatments.
A protocol for efficient induction of DNA damage under controlled in vitro conditions is not well established for barley.Our initial attempts to treat whole sterilized barley seeds according to established protocols in Arabidopsis failed.The reasons were high frequency of fungal contamination, highly variable responses to chemical treatments and complexity of root phenotype (multiple roots emerging).Based on extensive testing and modifications, we propose a protocol using dissected mature embryos and their cultivation on solid media containing genotoxins.Dissecting embryos represents an additional and demanding experimental step, but it helped in several ways.First, it greatly reduced fungal contaminations that occurred on whole seeds, possibly due to the complex surface of barley seeds.Second, it excluded endosperm which provided energy for the germinating embryo, allowing for minimal contact of the roots with the zeocin-containing media.Third, the germination showed less variation between the individuals and the experiments.
To induce DNA damage, we applied zeocin, a phleomycin D1/bleomycin-type antibiotic with radiomimetic effects.While DSBs are commonly caused by radiomimetic drugs, and bleomycin has been shown to induce DNA fragmentation 34 , research also revealed that only about 10% of DNA alterations caused by bleomycin involve both strands, and within this percentage, only a fraction are actual DSBs 35 .The ratio of SSBs to DSBs caused by bleomycin varies depending on concentration and chromatin compaction, ranging from 3:1 to 20:1 18,36,37 .Presence of both DSBs and SSBs was confirmed on DNA from bleomycin-treated barley root tissues via Comet assays 38 .Therefore, zeocin is a relatively broad-spectral DNA damage inducer that can potentially serve as a genotoxin for testing sensitivity of a wide range of DDR mutants.Induction of SSBs and to a lower extent also other non-DSB types of DNA damage by phleomycin type antibiotics 27 also justifies use of these chemicals for treatments of atr mutants.Although we did not assess the amount of DNA-SBs directly in our experiments, inhibitory effects of zeocin treatments on plant growth and root tissue differentiation indicated that it was effective already at the lowest applied dose (100 µg/ml).This was further supported by up-regulation of the DDR marker genes BRCA1 and RAD51.The effects on the shoot were less prominent, most likely because they were not in a direct contact with the drug-containing media.
Analysis of ATR mutant plants confirmed shorter but more seminal roots in hvatr.gcompared to WT 21 .This is a principal difference compared to phenotypic analysis of dicots, where the primary root offers a simple proxy for the plant growth under both mock and genotoxic stress conditions.The exact reason for this phenotype is unclear but we propose several speculative models to it.Presumed replication-coupled defects in hvatr.gmutant plants might shift the ration between the root elongation by active cell division versus initiation of the new roots.In wild-type plants, formation of too many new roots might be suppressed via interaction of ATR with the factors that integrate DDR with cell cycle such as SOG1 and WEE1 7,39,40 .Yet another mechanism might include metabolic problems because hvatr.g plants exhibited a down-regulation of transmembrane transport genes, most notably the nitrate and nitrogen-compound transporters.We also observed down-regulation of 13 EXPANSIN genes, www.nature.com/scientificreports/which could potentially impact cell growth and root number in hvatr.gplants.EXPANSINS were part of a root number quantitative trait locus in barley, although no specific gene was conclusively confirmed as the causative factor 41 .We also searched for misregulated genes known to be involved in barley root development.The only such candidate was down-regulated PME5 (PECTINESTERASE 11), which has previously been linked to root length regulation 42 .However, to answer this in an unbiased way, a complex forward-directed suppressor screen would have to be carried out in the hvatr.gbackground, seeking double mutants that restore wild-type-like phenotype.Genome-wide transcriptomics upon DNA damage proved useful in exploring plant DDR 7,22,43 .Therefore, it was interesting to observe transcriptional changes upon zeocin treatment in barley.A prominent effect was a response to oxidative stress.This is consistent with the notion that bleomycin type antibiotics, including zeocin, mediate production of superoxides and free oxygen radicals 44,45 .The upregulated enzymes included GSTUs and UGTs that use glutathione to detoxify reactive oxygen species.Glutathione is a buffer protein used against redox active molecules.In plants, some GSTUs have also strong antioxidative roles, with some classes of GSTU having peroxidase activity 46 .Furthermore, UGTs catalyze activation reactions of most secondary metabolites, including antioxidative molecules by addition of sugar moieties 47,48 .
Many zeocin-induced transcriptional changes included genes directly or indirectly associated with the cell cycle and cell division.The down-regulation of CYCs and CDKs indicated an inhibition of regular cell cycle progression.Presence of DNA damage generally leads to a halted cell cycle to gain time necessary for the repair 49 .A strong down-regulation of core histone genes and replication factors indicated reduced DNA synthesis.Many histone genes reach the peak of expression during S phase, when massive amounts of histones are needed for packaging of newly synthetized DNA.Interestingly, the down-regulation included also two DDR-associated H2A.X histone variant genes.Phosphorylated form of H2A.X (γH2A.X) marks the sites of active DNA damage repair 50 .Because H2A.X is present at specific genomic positions under non-DNA damaging conditions 51 , it likely that H2A.X is transcriptionally-regulated in the cell cycle stage-dependent manner.Our data suggest (indirectly) that the H2A.X response to DNA damage occurs at the level of post-translational modifications.An obvious consequence of the response at the protein rather than transcriptional level would be a faster reaction time.This could be important in order to pause cell cycle and start the repair before more damage occurs either directly by the mutagen or indirectly by continuation of regular cellular processes.The second prominent trend indicated a shift from mitotic cycling towards endoreduplication.Endoreduplication is a modified cell cycle where G2 phase is followed by another S-phase instead of mitosis 52 .Endoreduplication is part of a standard plant developmental program but can be alleviated by stress 39,53 .In Arabidopsis, the process is controlled by the KIP-RELATED PROTEINS (KRPs) and SIAMESE-RELATED proteins (SMRs), where the first promote mitosis while the latter endoreduplication 54 .We found zeocin-induced up-regulation of barley SMR3.In Arabidopsis, SMRs suppress mitosis by repressing A and B type CYCs and B type CDKs and number of these factors was downregulated by zeocin treatment in barley.The other SMR activity is to promote endoreduplication cycle by suppressing expression in the signaling cascade consisting of CYCDD, CDKA, RBR1, E2Fs and FBL17 and leading to KRPs 54 .In barley, we found transcriptionally down-regulated ERF factor DEL1 and FBL17.Barley plants show generally low levels of endoreduplication in somatic tissues under normal conditions but the frequency is higher in specialized tissues such as some endosperm developmental stages 20 .Interestingly, there was a significant increase in the frequency of endoreduplicated nuclei in the roots of zeocin-treated plants.This is in agreement with the molecular signatures in our RNA-seq analysis and suggests that barley performs adjustments to its cell cycle in responses to genotoxic stress via evolutionarily conserved SMR-dependent pathway.
To focus RNA-seq on DDR components, we created list of barley DDR genes using protein homologies with Arabidopsis candidates 22 .Through selection, filtering and manual curation steps, we identified 421 barley candidates.Some, like the barley BRCA1 homolog had lower similarity to Arabidopsis and were added manually.The list also included homologs of Arabidopsis DDR response master regulator SOG1, and its closely related genes ANAC044 and ANAC085 26,55 .We found five members of this family in barley.Comparisons with rice and maize suggested that SOG1 is a single copy gene in barley.In contrast, the other four family members clustered with one clade with the Arabidopsis ANAC044, ANAC085 and rice SGL.This suggested that these copies represent barley SGLs.Interestingly, SGLs were more transcriptionally responsive to zeocin-induced DNA damage than SOG1 in barley.This aligns with Arabidopsis, where SOG1 shows minimal transcriptional changes to DNA damage, but ANACs are up-regulated 7 , and indicates an evolutionarily conserved regulation of SOG1 at the protein level and SGLs at the transcriptional level during DDR.Other up-regulated genes included several DNA damage repair factors and positive regulators of HR including RAD51, BRCA1, RecQL3, NSE4.A, GMI1, or BCP1.
Transcriptomic data of hvatr.gshowed activated DDR even under the mock conditions, consistent with previous findings of increased DNA damage in non-treated hvatr.gplants 21 .Up-regulated factors indicated greater utilization of mismatch and NER pathways in ATR mutants.Notably, moderately up-regulated genes in hvatr.gincluded RAD23 and CDC48.RAD23 is connected to proteolysis helping with the cell-cycle progression and stress response 56 .In yeast, Rad23 recognizes a complex of Cdc48 and Ubiquitin ligase E4, facilitating proteasomal degradation.Plant CDC48 homologs have similar roles in protein degradation 57 .The enrichment of genes involved in ubiquitinilation, along with down-regulation of E2 Ub-conjugating enzymes and RGLG2 Ub-ligase E3 in hvatr.g,suggests increased protein turnover.Alternatively, elevated CDC48 expression may contribute to methionine metabolism for cell defense in the nicotineamine pathway or facilitate de-condensation of centromeric heterochromatin and/or activation of rDNA genes as observed in Arabidopsis.These pathways potentially have support from other upregulated genes.Notably, FBL17 was down-regulated in hvatr.g,and Arabidopsis FBL17 mutants were identified as SUPPRESSORS OF ATR 1 sensitivity phenotype 33 , indicating a transcriptional regulatory feedback loop between ATR and FBL17.
Besides of the obvious transcriptional differences between hvatr.g and wild-type, it has to be noted that vast majority of the transcriptional response between both genotypes remained unchanged.This is most likely due to a partial functional redundancy of plant ATM and ATR kinases as described in Arabidopsis 58  www.nature.com/scientificreports/barley atm atr double mutant would be needed to uncover the whole spectrum of genes controlled by these key DDR kinases.Alternatively, analysis of barley sog1 mutant might be methodologically easier and possibly even more informative option 7,8 In conclusion, our study indicates that barley exhibits a conserved response to chemically-induced DNA-SBs.We observed molecular signatures of oxidative stress response, that is consistent with the zeocin expected mode of action and responses to both DNA single and double strand breaks.We also identifies some genes that could be possible targets of modifying mitotic division and endoreduplication in barley.This study offers valuable resources for further detailed investigations into barley's DDR, its associations with other stresses, and plant development.

Plant materials and growth conditions
We used barley cultivars Golden Promise and Sebastian (WT), and hvatr.gTILLING mutant allele of ATR (HORVU.MOREX.r3.7HG0748510) gene 21 .In vitro plant cultivation was done in an air-conditioned phytochamber with a long day regime (16 h light, 150 µmol m −2 s −1 , 21 °C, 8 h dark, 19 °C).Plants used for seed production were grown in the climate-controlled phytotron under long day conditions (16 h light, 150 μmol m −2 s −1 , 15 °C, 8 h dark, 13 °C, 65% humidity).Plant materials used in this work are a part of cultivated cereal varieties and were not resourced in field.Golden Promise (Acc.No. HOR 16645) is available at the Leibniz Institute of Plant Genetics and Crop Plant Research, IPK, Genebank, Gatersleben, Germany and Sebastian (Acc.No. 03C0602773) at the Germplasm Resource Information Network (GRIN), Prague, Czech Republic.Mutant hvatr.g is available at HorTILLUS (Hordeum-TILLING-University of Silesia) database upon request from M.S. Experimental research on plant material in this study, including its collection, complied with the relevant institutional, national, and international guidelines and legislation.

DNA damage assays
Barley seeds were surface sterilized with 70% ethanol (v/v) for 5 min, followed by 10 min treatment with 8% sodium hypochlorite (v/v) and final triple wash with ddH 2 O. Sterilized seeds were imbibed overnight in sterile water at 4 °C in dark.The following day, embryos were carefully excised under binocular in a sterile laminal flow-hood from the seeds and placed scutellum side down on ½ MS medium with 0.6% agarose (w/v) or medium supplemented with zeocin (ThermoFisher Scientific, Cat.no.R25001) in 107 × 94 × 96 mm boxes (Duchefa, Cat.no.S1686).In place of a lid, another container was sealed to the one containing medium with embryos with a parafilm (total height 192 mm).Containers were placed in a phytochamber (Percival Scientific) and plants were grown for 14 days after which they were carefully pulled out of the medium for measurements.All the measurements were completed using the ImageJ plugin SmartRoot 59 .Experiments were performed in three biological replicates.For testing the different concentrations of zeocin required for wild-type barley (cv.Golden Promise) phenotypic response each replicate contained 18-20 individual plants.In the DNA damage assays comparing the phenotypic response to zeocein of hvatr.g and wild-type cv.Sebastian, each replicate had 9-10 plants per genotype.Statistical significance was tested with One-way ANOVA with post hoc Tukey HSD in Minitab (www.minit ab.com).

Propidium iodide staining and root microscopy
Barley plants cv.Golden Promise grown for seven days on a control or 100 μg/ml zeocin media were used to assess root morphology at microscopic level.Three representative plant samples were chosen for imaging.The whole root was stained by pseudo-Schiff propidium iodide staining as described (Coiro and Truernit 60 ).Incubation with propidium iodide and Schiff reagent lasted for 48 h.Following the overnight incubation with chloral hydrate solution roots were mounted on glass in water.Imaging was performed using Leica TCS SP8 STED3X confocal microscope (Leica Microsystems, Wetzlar, Germany), equipped with an HC PL APO CS2 10×/0,4 DRY objective, hybrid detectors (HyD), and the Leica Application Suite X (LAS-X) software version 3.5.5 with the Leica Lightning module (Leica, Buffalo Grove, IL, USA).For propidium iodide acquisition, laser excitation at 534 nm and emission at 550-730 nm were used.The maximal projection pictures were constructed from aligned Z-stack images of approximately 250-300 μm steps, containing 45 individual optical sections.The images were post-processed by Leica Lightening software module.

Flow-cytometry
The nuclear ploidy measurements were done on 20 Golden Promise plants grown from dissected embryos for seven days on solid ½ MS medium with or without 300 μg/ml zeocin.Two to three root apical meristems from each individual plant were chopped using a razor blade directly into 500 ml Otto I buffer (0.1 M citric acid, 0.5% Tween 20 v/v).Nuclei suspension was filtered through a 50 µm nylon mesh into a fresh tube, mixed with 1 ml of Otto II buffer (0.

Reverse transcription-quantitative polymerase chain reaction (RT-qPCR)
Golden Promise seeds were surface sterilized, embryos dissected, and grown on ½ MS medium with 0.6% agarose.After 7 days, young seedlings were moved to liquid ½ MS medium without or with zeocin (300 μg/ ml).Sampling was performed at strict time points 0.5, 1 and 6 h after the beginning of treatment.All root apical meristems (RAMs) from a single plant were dissected and flash-frozen in liquid nitrogen.Three plants were taken for each treatment and stored until use at − 80 °C.RNA extraction was performed by RNeasy Mini Kit

Figure 1 .
Figure 1.Wild-type barley (cv.Golden Promise) parameters in response to zeocin treatment.(A) Representative phenotypes of barley seedlings grown for 7 days on ½ MS (mock) and zeocin concentrations (100, 300, 500 μg/ml) containing media.Scale bars = 5 cm.(B) Effect of zeocin treatments on barley seedling growth.Error bars indicate the standard deviation between the means of three biological replicates.The letters above columns indicate similarities between samples.The same letters indicate samples that were not significantly different in one-way ANOVA with post hoc Tukey's test (P < 0.05).Kruskall-Wallis H-test with post hoc Conover-Iman test of multiple comparisons with Benjamini-Hochberg procedure (P < ½ α, α = 0.05) were used to analyze differences in root numbers.NS, not significant.(C) Representative confocal microscopy maximal projection images of the root surface in plants grown in mock conditions and plants treated with 100 μg/ml zeocin.Scale bar = 200 μm.(D) Percentage of nuclei at specific cell cycle stages found in the root tips of plants grown on mock and medium with 300 μg/ml zeocin.Super cycle value describing the amount of nuclei in endoreduplication (END).Statistical relevance distinguished by two-sample T-test, (P < 0.05).*P < 0.05, ***P < 0.001.

Figure 2 .
Figure 2. Barley atr mutant plants are sensitive to zeocin treatment.(A) Representative images of wild-type cultivar Sebastian (Seb) and atr mutant plants treated with 100 μg/ml zeocin, or untreated (mock).Scale bars = 5 cm.(B) Effect of zeocin treatment on Seb and atr plants' root growth.Error bars indicate the standard deviation between the means of three biological replicates.For the absolute values represented the letters above columns indicate similarities between samples.The same letters indicate samples that were not significantly different in one-way ANOVA with post hoc Tukey's test (P < 0.05).For the assessment of the statistical significance of the relative values the two-sample T-test was used, *P < 0.05, **P < 0.01.(C) Number of roots in wild-type (Seb) and atr plants grown in mock and genotoxic conditions.Kruskall-Wallis H-test with post hoc Conover-Iman test of multiple comparisons with Benjamini-Hochberg procedure (P < ½ α, α = 0.05) was used to analyze differences in root numbers.NS, not significant, ***P < 0.001.(D) Effect of zeocin treatment on Seb and atr plants' shoot growth.Error bars indicate the standard deviation between the means of three biological replicates.For the absolute values represented the letters above columns indicate similarities between samples.The same letters indicate samples that were not significantly different in one-way ANOVA with post hoc Tukey's test (P < 0.05).For the assessment of the statistical significance of the relative values the two-sample T-test was used, *P < 0.05.

Figure 3 .
Figure 3. Phylogenetic analysis of barley SOG1 and SGLs.(A) Multiple sequence alignment of NAC protein domains present in five barley SOG1 protein candidates with the NAC domains of Arabidopsis thaliana (At) SOG1, ANAC044, and ANAC085, Oryza sativa (Os) SOG1 and SOG1-like (SGL), and Zea mays (Zm) NACTF8 and NACTF99 proteins by MUSCLE.Amino acids are highlighted based on consensus sequence and their physicochemical properties.(B) Phylogenetic Maximum Likelihood tree built based of the multiple sequence alignment, with Arabiodopsis ANAC005 and ANAC006 NAC protein domains used as an outgroup.Bootstrap values are shown next to the branches, distance scale = 0.5.

Figure 4 .
Figure 4. Barley BRCA1 and RAD51 are transcriptionally upregulated after zeocin treatment in Golden Promise wild-type.Reverse transcription-quantitative PCR values were normalized to the expression of UBC2 as mean normalized expression (MNE).Error bars indicate standard deviation of three biological replicates.Statistical significance of the difference presented was tested by two-sample T-test, (P < 0.05).*P < 0.05.

Figure 5 .
Figure5.Zeocin effect on gene expression in wild-type barley.(A) Biological processes enriched among significantly up-regulated and down-regulated genes in zeocin versus mock treated wild-type (Sebastian) plants.Redundant GO terms were removed manually, based on P-value.The full list of GO terms can be found in Suppl.Dataset 3. Statistical significance was determined by Fisher's one-tailed test with g:SCS algorithm correction.Gene ratio represents the number of genes found in the category compared to the total number of genes in the query.(B) Volcano plot representing all genes detected as differentially expressed in wild-type (Sebastian) plants treated with zeocin relative to mock treatment.Horizontal dashed line indicates genes passing P < 0.05.Vertical dashed lines separate genes with log 2 Fold Change ≤ − 1 (blue) or log 2 Fold Change ≥ 1 (red).Specific marker genes were highlighted by names.

Figure 6 .
Figure 6.Effects of ATR mutation on gene expression in barley.(A) Gene enrichment analysis for biological processes in atr mutant compared to wild-type (Sebastian).Redundant GO terms were removed manually, based on P-value.The full list of GO terms can be found in Suppl.Dataset 3. Statistical significance was determined by Fisher's one-tailed test with g:SCS algorithm correction.Gene ratio represents the number of genes found in the category compared to the total number of genes in the query.(B) Volcano plot representing all genes detected as differentially expressed in atr plants relative wild-type (Sebastian) plants.Horizontal dashed line indicates genes passing P < 0.05.Vertical dashed lines separate genes with log 2 Fold Change ≤ − 1 (blue) or log 2 Fold Change ≥ 1 (red).Specific marker genes were highlighted by names.

Figure 7 .
Figure 7. Response of atr mutant plants to zeocin treatment.(A) Volcano plot representing all genes differentially expressed in zeocin-treated versus mock-treated atr mutant plants.Horizontal dashed line indicates genes passing P < 0.05.Vertical dashed lines separate genes with log 2 Fold Change ≤ − 1 (blue) or log 2 Fold Change ≥ 1 (red).Specific marker genes were highlighted by names.(B) Venn diagrams showing differentially expressed genes after zeocin treatment in wild-type and atr mutant plants relative to the mock treatment of the same genotype. https://doi.org/10.1038/s41598-024-53264-0